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1 Introduction 

It has been noticed in recent years that a significant non-trivial physical phenomena 
in granular gases can be described mathematically by dissipative Boltzmann type 
equations, as can be seen in [17] for a review in the area. As motivated by this 
particular phenomena of energy dissipation at the kinetic level, we consider in this 
chapter the Boltzmann equation for non-linear interactions of Maxwell type and 
some generalizations of such models. 

The classical conservative (elastic) Boltzmann equation with the Maxwell-type 
interactions is well-studied in the literature (see [5, 14] and references therein). 
Roughly speaking, this is a mathematical model of a rarefied gas with binary colli- 
sions such that the collision frequency is independent of the velocities of colliding 
particles, and even though the intermolecular potentials are not of those correspond- 
ing to hard sphere interactions, still these models provide a very rich inside to the 
understanding of kinetic evolution of gases. 
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Recently, Boltzmann equations of Maxwell type were introduced for models of 
granular gases were introduced in [7] in three dimensions, and a bit earlier in [3] for 
in one dimension case. Soon after that, these models became very popular among 
the community studying granular gases (see, for example, the book [13] and ref- 
erences therein). There are two obvious reasons for such studies The first one is 
that the inelastic Maxwell-Boltzmann equation can be essentially simplified by the 
Fourier transform similarly as done for the elastic case, where its study becomes 
more transparent [6,7]. The second reason is motivated by the special phenomena 
associated with homogeneous cooling behavior, i.e. solutions to the spatially ho- 
mogeneous inelastic Maxwell-Boltzmann equation have a non-trivial self-similar 
asymptotics, and in addition, the corresponding self-similar solution has a power- 
like tail for large velocities. The latter property was conjectured in [16] and later 
proved in [9, 11]. This is a rather surprising fact, since the Boltzmannn equation 
for hard spheres inelastic interactions has been shown to have self similar solutions 
with all moments bounded and large energy tails decaying exponentially. The con- 
jecture of self-similar (or homogeneous cooling) states for such model of Maxwell 
type interactions was initially based on an exact Id solution constructed in [1]. It is 
remarkable that such an asymptotics is absent in the elastic case (as the elastic Boltz- 
mann equation has too many conservation laws). Later, the self-similar asymptotics 
was proved in the elastic case for initial data with infinite energy [8] by using an- 
other mathematical tools compared to [9] and [12]. 

Surprisingly, the recently published exact self-similar solutions [12] for elastic 
Maxwell type model for a slow down process, derived as a formal asymptotic limit 
of a mixture, also is shown to have power-like tails. This fact definitely suggests 
that self-similar asymptotics are related to total energy dissipation rather than local 
dissipative interactions. As an illustration to this fact, we mention some recent pub- 
lications [2, 15], where Id Maxwell-type models were introduced for non-standard 
applications such as models in economics and social interactions, where also self- 
similar asymptotics and power-like tail asymptotic states were found. 

Thus all the above discussed models describe qualitatively different processes 
in physics or even in economics, however their solutions have a lot in common 
from mathematical point of view. It is also clear that some further generalizations 
are possible: one can, for example, include in the model multiple (not just binary) 
interactions still assuming the constant (Maxwell-type) rate of interactions. Will 
the multi-linear models have similar properties? The answer is yes, as we shall see 
below. 

Thus, it becomes clear that there must be some general mathematical properties 
of Maxwell models, which, in turn, can explain properties of any particular model. 
That is to say there must be just one main theorem, from which one can deduce 
all above discussed facts and their possible generalizations. Our goal is to consider 
Maxwell models from very general point of view and to establish their key proper- 
ties that lead to the self-similar asymptotics. 

All the results presented in this chapter are mathematically rigorous. Their full 
proofs can be found in [10]. 
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After this introduction, we introduce in Section [2] three specific Maxwell mod- 
els of the Boltzmann equation: (A) classical (elastic) Boltzmann equation; (B) the 
model (A) in the presence of thermostat; (C) inelastic Boltzmann equation for 
Maxwell type interactions. Then, in Section [3] we perform the Fourier transform 
and introduce an equation that includes all the three models as particular cases. A 
further generalization is done in Section @] where the concept of generalized multi- 
linear Maxwell model (in the Fourier space) is introduced. Such models and their 
generalizations are studied in detail in Sections[5]and [6] The most important for our 
approach concept of L-Lipschitz nonlinear operator is explained in Section [4] It is 
shown (Theorem 4.2) that all multi-linear Maxwell models satisfy the L-Lipschitz 
condition. This property of the models constitutes a basis for the general theory. 

The existence and uniqueness of solutions to the initial value problem is stated 
in Subsection l5. 1 K Theorem 5.2). Then, in Subsection l5.2l we present and study the 
large time asymptotics under very general conditions that are fulfilled, in particu- 
lar, for all our models. It is shown that L-Lipschitz condition leads to self-similar 
asymptotics, provided the corresponding self-similar solution does exist. The exis- 
tence and uniqueness of such self-similar solutions is stated in Subsection l5 .3 K The- 
orem 5.12). This theorem can be considered, to some extent, as the main theorem 
for general Maxwell-type models. Then, in Subsection 15.41 we go back to multi- 
linear models of Section [4] and study more specific properties of their self-similar 
solutions. 

We explain in Section [6] how to use our theory for applications to any specific 
model: it is shown that the results can be expressed in terms of just one function 
H (p) , p > 0, that depends on spectral properties of the specific model. General prop- 
erties (positivity, power-like tails, etc.) self-similar solutions are studied in Subsec- 
tions [6j]and[62] It includes also the case of Id models, where the Laplace (instead 
of Fourier) transform is used. In Subsection 6.3, we formulate, in the unified state- 
ment (Theorem 11.1), the main properties of Maxwell models (A),(B) and (C) of 
the Boltzmann equation. This result is, in particular, an essential improvement of 
earlier results of [7] for the model (A) and quite new for the model (B). 

Applications to one dimensional models are also briefly discussed at the end of 
Subsection 6.3. 



2 Maxwell models of the Boltzmann equation 

We consider a spatially homogeneous rarefied cZ-dimensional gas (d = 2,3, . . .) of 
particles having a unit mass. Let /(v, t ), where v G W 1 and t 6 R + denote respectively 
the velocity and time variables, be a one-particle distribution function with usual 
normalization 



Then f(v,t) has an obvious meaning of a time-dependent probability density in 
M. d . We assume that the collision frequency is independent of the velocities of the 




(1) 



4 A.V. Bobylev, C. Cercignani and I.M. Gamba 

colliding particles (Maxwell-type interactions). We discuss three different physical 
models (A), (B) and (C). 

(A) Classical Maxwell gas (elastic collisions). In this case f(v,t) satisfies the 
usual Boltzmann equation 

ft = QW,f) = /„ dwda>g(^-)\f{v')f{w')-f{v)f{w)] , (2) 
Jw'xs d -' \u\ 

where the exchange of velocities after a collision are given by 

v' = -(v + w+ \u\co) , and w' = —(v + w— \u\co) 

where u = v — w is the relative velocity and Q € S . For the sake of brevity we 
shall consider below the model non-negative collision kernels g(s) such that g(s) is 
integrable on [—1,1]. The argument t of f(v,t) and similar functions is often omitted 
below (as in Eq. (01). 

(B) Elastic model with a thermostat This case corresponds to model (A) in the 
presence of a thermostat that consists of Maxwell particles with mass m > having 
the Maxwellian distribution 

Mv= )-"/2 exp( __LL) (3 ) 

m 11 

with a constant temperature T > 0. Then the evolution equation for/(^,f) becomes 

ft = Q(f,f) + [ dwd(Og{^-)\f{v')M{w')-f{v)M{w)] , (4) 
J \u\ 

where 6 > is a coupling constant, and the exchange of velocities is now 

, v + m(w + \u\(o) , , v+mw—\u\(0 

v = — , and w = — , 

1 +m l+m 

with it = v — w the relative velocity and , (0 E 

Equation (01 was derived in [12] as a certain limiting case of a binary mixture of 
weakly interacting Maxwell gases. 

(C) Maxwell model for inelastic particles. We consider this model in the form 
given in [9]. Then the inelastic Boltzmann equation in the weak form reads 



-(f, V )= [ dvdwdcof(v)f(w) ^1 [y/(v') - V (y)] , (5) 

at j VL d xR d xS d-\ \ u \ 

where y/(v) is a bounded and continuous test function, 

if,Y)= [dvf(v,t)y(v),u = v-w,(0£S d ~ 1 ,v' = v-^4^-{u-co)co, (6) 
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the constant parameter < e < 1 denotes the restitution coefficient. Note that the 
model (C with e = 1 is equivalent to the model (A) with some kernel g(s). 

All three models can be simplified (in the mathematical sense) by taking the 
Fourier transform. 

We denote 

f(k,t) = -nf] = (f,e- ik - v ), k e R d , (7) 

and obtain (by using the same trick as in [6] for the model (A)) for all three models 
the following equations: 

(A) ft=Q(f,f) = J s dco g { k -^-)[f{k + )f{k_)-f{k)fm , 

(8) 

where k± = \(k± \k\co), (0 e S*' 1 , /(0) = 1. 

(B) f t = Q(fJ) + eJ s dag(^){f(k + )M(k-)-f(k)M(0)} , 

(9) 

where M(jfc) = e~ # , k+ = k -±^ , k- = k-k+, (0 6 S d ~\ /(0) = 1 . 

(C) /, = J s dJ-^\f(k + )f(k-)-f(k)fm , 

(10) 

where /(0) = 1, k+ = ^(k ■ G))(D, k-=k— k+, with co G S d ~ l is the direction 
containing the two centers of the particles at the time of the interaction. Equivalently, 
one may alternative write k- = ^-(k — \k\na>), and k + = k — k-, where now (b € 

S d ~ l is the direction of the post collisional relative velocity, and the term dw is 
replaced by a function g(^-)da). 

Case (B) can be simplified by the substitution 

~ T\k\ 2 
f(k,t)=f(k,t)exp{ LL] , (11) 

leading, omitting tildes, to the equation 

(B') f t =Q(f,f) + e Ucog( k ^)[f( k -±^)-nk)] , 

Js*- 1 \k\ l+m 

(12) 

i.e., the model for (B) with T = 0, or equivalently a linear collisional term the back- 
ground singular distribution. Therefore, we shall consider below just the case (B'), 
assuming nevertheless that f(k, t) in Eq. ( fT2l is the Fourier transform (0 of a prob- 
ability density f(v,t). 
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3 Isotropic Maxwell model in the Fourier representation 



We shall see that these three models (A), (B) and (C) admit a class of isotropic so- 
lutions with distribution functions / = /(|v| , t ). Indeed, according to ^ we look for 
solutions / = f(\k\,t) to the corresponding isotropic Fourier transformed problem, 
given by 

x=\k\ 2 , (p(x,t)=f(\k\,t)=&\f(\v\,t)}, (13) 
where (p(x,t) solves the following initial value problem 



(p,= dsG(s){(p[a(s)x](p[b(s)x]-(p(x)}- 



d S H( S ){(p[c( S )x}-cp(x)} 

<Pt=o = <poW, <p(o,0 = i ; 



(14) 



where a(s), b(s), c(s) are non-negative continuous functions on [0, 1], whereas G(s) 
and H(s) are generalized non-negative functions such that 



dsG(s) < °o 



dsH(s) < oo 



(15) 



Thus, we do not exclude such functions as G = S(s — so), < sq < 1, etc. We shall 
see below that, for isotropic solutions (1131) . each of the three equations (0, dXCTb . ( fT~2T > 
is a particular case of Eq. (114b , 

Let us first consider Eq. (© with f(k,t) = (p(x,t ) in the notation (1131 1. In that case 



and the integral in Eq. $Q reads 

d(Qg((O - (0)(p 



£0o = — eS' 1 - 1 ,d = 2, 



\k\ 



sp- 



it is easy to verify the identity 



1 + COq-CO 



l-coo-co 



(16) 



d(0F((0 ■ coo) = \S d ~ 2 \ I dzF(z)(l-z 



S d-i 



(17) 



where l^ -2 ! denotes the "area" of the unit sphere in R d ~ l for d > 3 and | ^S" 1 = 2. 
The identity ( fTTI ) holds for any function F(z) provided the integral as defined in the 
right hand side of dTVb exists. 
The integral (fTol l now reads 
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|S"- 2 | |Y^( Z )(1- Z 2)^ ( p(xl±^)(p(xi y i 







dsG(s)<p(sx)<p[(l - 



six 



where 

G(s)=2 d - 2 \S d - 2 \g(l-2 S )[s(l-s)}^ , d = 2,3,.... (18) 
Hence, in this case we obtain Eq. ( TBI ), where 

(A) a(s)=s, b(s) = l-s, H{s)=0, 

(19) 

G(s) is given in Eq. (fT8l ). 

Two other models (B') and (C), described by eqs. (1121) . (TToT > respectively, can be 
considered quite similarly. In both cases we obtain Eq. ( TBI , where 

4/72 

(B') a(s) = s, b(s) = 1-5, c(s) = 1 - — — -js, 

(1 + ;«J Z 

#0) = 0G(s) , 

(20) 

G(s) is given in Eq. (fT8l : 

(O «W = ^, *(*) = ! (1+e)(3 - e) " 



4 ' w 4 



H(*)=0, G(.) = ^'- 2 |(l 



rf-3 



(21) 



Hence, all three models are described by Eq. ( TBI ) where < a(s),b(s),c(s) < 1 
are non-negative linear functions. One can also find in recent publications some 
other useful equations that can be reduced after Fourier or Laplace transformations 
to Eq. (TBI) (see, for example, [2], [15] that correspond to the case G = 8(s — so), 
H = 0). 

The equation (TBI with H (s) = first appeared in its general form in [9] in con- 
nection with models (A) and (C). The consideration of the problem of self-similar 
asymptotics for Eq. ( fT4b in that paper made it quite clear that the most important 
properties of "physical" solutions depend very weakly on the specific functions 
G(s), a(s) and b(s). 



4 Models with multiple interactions 



We present now a general framework to study solutions to the type of problems 
introduced in the previous section. 
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We assume, without loss of generality, (scaling transformations t = at, a = 
const.) that 

f ds[G(s)+H(s)] = l (22) 
Jo 

in Eq. ( TT4b . Then Eq. (fl4] i can be considered as a particular case of the following 
equation for a function u{x,t) 

u t + u = r(u), x>0, f>0, (23) 

where 



r(u) = £ a„r^{u) , £ a„ = 1 , a,, > , 

n=\ n=\ ^ (24) 

r i - n \u) = jdai... jda„A„(ai,...,a„)Y[u(a k x), n = l,...,N. 
We assume that 

A„(a) =A n (ai,...,a n ) > , / da\...\ da„ A(a\, . . . ,a n ) = 1 , (25) 

Jo Jo 

where A„(a) = A„(ai , . . . ,a„) is a generalized density of a probability measure in 
for any n = 1 , . . . ,7Y. We also assume that all A„ (a) have a compact support, i.e., 

n 

A n ( ai ,...,a„) = if Y, a l > r2 > n = l,...,N, (26) 

k=l 

for sufficiently large < R < °°. 

Eq. (O is a particular case of Eq. d23l with 



N = 2 , a\ = I dsH(s) , c«2 



1 



[ l dsG(s) 
Jo 



Ai(ai) = — / rfi//(j)5[ai -c(s)l (27) 
a\ Jo 

1 r 1 

A 2 (ai,a2) = — / dsG(s)8[a\-a(s)]5[a 2 -b(s)] . 
0C2 Jo 

It is clear that Eq. ( |23l can be considered as a generalized Fourier transformed 
isotropic Maxwell model with multiple interactions provided u(Q,t) = 1, the case 
N = °° in Eqs. (|24| | can be treated in the same way. 
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The general problem we consider below can be formulated in the following way. 
We study the initial value problem 

u t + u = r(u) , ki (= o = wo(x) , x>0, t>0, (28) 

in the Banach space B = C(M+) of continuous functions u(x) with the norm 

||w|| = sup \u(x)\ . (29) 

x>0 

It is usually assumed that ||uo|| < 1 and that the operator F is given by Eqs. 
d24"l) . On the other hand, there are just a few properties of F(m) that are essential 
for existence, uniqueness and large time asymptotics of the solution u(x,t) of the 
problem d28l . Therefore, in many cases the results can be applied to more general 
classes of operators F in Eqs. ( 128b and more general functional space, for example 
B = C(W l ) (anisotropic models). That is why we study below the class (l24l of oper- 
ators F as the most important example, but simultaneously indicate which properties 
of F are relevant in each case. In particular, most of the results of Section 4-6 do 
not use a specific form (l24l i of F and, in fact, are valid for a more general class of 
operators. 

Following this way of study, we first consider the problem (f28j) with F given by 
Eqs. (1241) and point out the most important properties of F. 

We simplify notations and omit in most of the cases below the argument x of the 
function u(x,t). The notation u(t) (instead of u(x,t)) means then the function of the 
real variable t > with values in the space B = C(K+). 

Remark 1. We shall omit below the argument x £ M+ of functions u(x), v(x), etc., 
in all cases when this does not cause a misunderstanding. In particular, inequalities 
of the kind \u\ < |v| should be understood as a point- wise control in absolute value, 
i.e. "|m(x)| < |v(x) for any x > 0" and so on. 

We first start by giving the following general definition for operators acting on a 
unit ball of a Banach space B denoted by 

U = {u eB: \\u\\ < 1} (30) 

Definition 1. The operator F = F (u) is called an L-Lipschitz operator if there exists 
a linear bounded operator L.B^B such that the inequality 

|r(«i)-r(«2)|<L(|wi-H2|) (3i) 

holds for any pair of functions M12 in U. 

Remark 2. Note that the L-Lipschitz condition PTT ) holds, by definition, at any point 
x £ (or x E W 1 if B = C(M. d )). Thus, condition OTT l is much stronger than the 
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classical Lipschitz condition 

||r(ui)-r(u 2 )||<c||iii-ii2|| if u l2 eu (32) 

which obviously follows from d3~TT i with the constant C = \\L\\b, the norm of the 
operator L in the space of bounded operators acting in B. In other words, the ter- 
minology "L-Lipschitz condition" means the point-wise Lipschitz condition with 
respect to an specific linear operator L. 

We assume, without loss of generality, that the kernels A n [a y, . . . ,a„) in Eqs. d24i i 
are symmetric with respect to any permutation of the arguments (a\,...,a n ), n = 
2,3,. ..,N. 

The next lemma states that the operator r(u) defined in Eqs. ([24-b . which satis- 
fies r(l) = 1 (mass conservation) and maps U into itself, satisfies an L-Lipschitz 
condition, where the linear operator L is the one given by the linearization of F near 
the unity. See [10] for its proof. 

Theorem 1. The operator T (u) defined in Eqs. d24b maps U into itself and satisfies 
the L-Lipschitz condition ( 1311 ), where the linear operator L is given by 



Lu = J daK(a)u(ax) , (33) 

with 



(34) 



N 

K ( a ) = Y, na n K "( a )> 
;i=l 

,00 N 

where K„{a) = I dai--- \ da„ A„(a,a2, ■ ■ ■ ,a„) and } CC n = \ . 

Jo Jo , I=1 

for symmetric kernels A n (a,a2, ...,a n ),n = 2, ...N . 

And the following corollary holds. 

Corollary 1. The Lipschitz condition ( 1321 ) is fulfilled for r(u) given in Eqs. ( 1241 ) 
with the constant 

N 

C=||L|| = £«a„, £a„ = l, (35) 

n=] n=l 

where \\L\\ is the norm ofL in B. 

It can also be shown that the L-Lipschitz condition holds in B = C(R ) for "gain- 
operators" in Fourier transformed Boltzmann equations (0, (0 and ( fTOb . 
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5.1 Existence and uniqueness of solutions 

It is possible to show, with minimal requirements, the existence and uniqueness 
results associated with the initial value problem d28l ) in the Banach space (B, \\ ■ 
||), where the norm associated to B is defined in ( |29l l In fact, this existence and 
uniqueness result is an application of the classical Picard iteration scheme and holds 
for any operator F which satisfies the usual Lipschitz condition (f32t and transforms 
the unit ball U into itself. The proof of all statements below can be found in [10]. 

Lemma 1. (Picard Iteration scheme) The operator T (u) maps U into itself and sat- 
isfies the L-Lipschitz condition ( 1321 ), then the initial value problem ( 1281 l with arbi- 
trary uq GU has a unique solution u(t) such that u{t) £ U for any t > 0. 

Next, we observe that the L-Lipschitz condition yields an estimate for the differ- 
ence of any two solutions of the problems in terms of their initial states. This is a 
key fact in the development of the further studies of the self similar asymptotics. 

Theorem 2. Consider the Cauchy problem ( 128b with \\uq\\ < 1 and assume that the 
operator T : B — > B 

(a) maps the closed unit ball U C B to itself and 

(b) satisfies a L-Lipschitz condition QTl ) for some positive bounded linear operator L : 
B^B. 

Then 

i) there exists a unique solution u{t) of the problem ( 1281 ) such that ||m(?)|| < I for any 
t > 0; 

ii) any two solutions u{t) and w{t) of problem d28l ) with initial data in the unit ball U 
satisfy the inequality 

\u(t)-w{t)\<exp{t(L-l)}(\u -wo\) . (36) 

Note that under the same conditions as in Theorem [T] the operator T given in 
Eqs. d24b satisfies necessary conditions for the Theorem|2] 

We remind to the reader that the initial value problem (l28l i appeared as a general- 
ization of the initial value problem ( TBI ) for a characteristic function (p(x,t ), i.e., for 
the Fourier transform of a probability measure (see Eqs. ( fT~3b , (Q])). It is important 
therefore to show that the solution u(x, t) of the problem d28l ) is a characteristic func- 
tion for any t > provided this is so for t = 0, which is addressed in the following 
statement. 
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Lemma 2. Let U' C U C B be any closed convex subset of the unit ball U (i.e., 
u = (1 - d)u\ + 9u 2 G U' for any mi, 2 G {/' and 9 G [0, 1]). //"mo G £/' in Eg. (f28b 
and/ {/ is replaced by U' in the condition (1) of Theorem [2] ?/ie theorem holds and 
u(t) G U' for any t > 0. 

Remark 3. It is well-known (see, for example, the textbook [18]) that the set U' C U 
of Fourier transforms of probability measures in Mr (Laplace transforms in the case 
of M_|_) is convex and closed with respect to uniform convergence. On the other hand, 
it is easy to verify that the inclusion F(I/') C U', where F is given in Eqs. (l24l) . 
holds in both cases of Fourier and Laplace transforms. Hence, all results obtained 
for Eqs. d23l ). d24l) can be interpreted in terms of "physical" (positive and satisfying 
the condition (fl}) solutions of corresponding Boltzmann-like equations with multi- 
linear structure of any order. 

We also point out that all results of this section remain valid for operators F 
satisfying conditions (T24l) . with a more general condition such as 

N 

£ On < 1 , On > , (37) 

n=l 

so that F(l) < 1 and so the mass may not be conserved. The only difference in this 
case is that the operator L satisfying conditions d3~3l >. d34l > is not a linearization of 
r(u) near the unity, but nevertheless Theorem[T]remains true. The inequality (|37| | is 
typical for Fourier (Laplace) transformed Smoluchowski-type equations where the 
total number of particles is decreasing in time (see [21, 22] for related work). 



In the next three sections we study in more detail the solutions to the initial 
value problem (l28li-(|2"9ii constructed in Theorem[2]and, in particular, their long time 
behavior, existence, uniqueness, and properties of the self-similar solutions. 



5.2 Large time asymptotics 

The long time asymptotics results are a consequence of some very general properties 
of operators F, namely, that F maps the unit ball U of the Banach space B = C(R+) 
into itself, F is an L-Lipschitz operator (i.e. satisfies OU) and that F is invariant 
under dilations. 

These three properties are sufficient to study self-similar solutions and large time 
asymptotic behavior for the solution to the Cauchy problem d28T > in the unit ball U 
of the Banach space C(R+). 
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Main properties of the operator T: 

a) r maps the unit ball U of the Banach space B = C(K+) into itself, that is 

\\r(u) || < 1 for any «6C(R+) such that \\u\\<\. (38) 

b) r is an L-Lipschitz operator (i.e. satisfies OTb ) with L from d33l . i.e. 

|F(mi)— r{u2)\{x) <L{\u\ — M2DM = / daK(a)\u\(ax) —U2(ax)\ , 

Jo 

(39) 

for K(a) > 0, for all x > and for any two functions U\ 2 £ C(E + ) such that 

IK 2 ||<i. 

c) F is invariant under dilations: 

e x9 r{u)=r{e^u) , ^=*4-, e T ®«(x)= M (xe T ), Tel. (40) 

No specific information about F beyond these three conditions will be used in 
this section. 

It was already shown in Theorem [2] that the conditions a) and b) guarantee ex- 
istence and uniqueness of the solution u(x,t) to the initial value problem (l28l-d29]i. 
The property b) yields the estimate (l36l l that is very important for large time asymp- 
totics, as we shall see below. The property c) suggests a special class of self-similar 
solutions to Eq. ( f28l >. 

We recall the usual meaning of the notation y = 0(x p ) (often used below): y = 
0(x p ) if and only if there exists a positive constant C such that 

\y(x)\<Cx p for any x > 0. (41) 

In order to study long time stability properties to solutions whose initial data 
differs in terms of 0(x p ), we will need some spectral properties of the linear operator 
L. 

Definition 2. Let L be the positive linear operator given in Eqs. (1331 . (l34l.then 

Lx p = X(p)x p , 0<X(p)= [ daK(a)a p <°°, p>0, (42) 

Jo 

and the spectral function ju(p) is defined by 

M(P) = ^- (43) 

An immediate consequence of properties a) and b), as stated in (|39l . is that one 
can obtain a criterion for a point-wise in x estimate of the difference of two solu- 
tions to the initial value problem (l28l yielding decay properties depending on the 
spectrum of L, as the following statement and its corollary assert. 
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Lemma 3. Let « 1.2 be any two classical solutions of the problem ([28} with 
initial data satisfying the conditions 



|ki i2 (jc,0)| < 1, |ki(x,0)-h 2 (x,0)| <Cx'\ x>0 
for some positive constant C and p. Then 

\m(x,t)-u z (x,t)\ <Cx p e-' {l - l[p)) , forall t>0 
Corollary 2. The minimal constant C for which condition (144b is satisfied is 



C 



sup - 

x>0 



|ki( x,0) -u 2 (x,0)\ 



u\{x,0) —U2(x,0)\ 



xP 



and the following estimate holds 
ui(x,t) — U2(x,t)\ 



xP 



<e 



-t(i-l(p)) 



u\(x,0) —U2(x,0)\ 



xP 



(44) 



(45) 



(46) 



(47) 



for any p > 0. 



A result similar to Lemma[3]was first obtained in [9] for the inelastic Boltzmann 
equation whose Fourier transform is given in example (C), Eq. (TTOb . Its corollary 
in the form similar to (l47l i for equation ( TTOb was stated later in [11] and was inter- 
preted there as "the contraction property of the Boltzmann operator" (note that the 
left hand side of Eq.d47li can be understood as a distance between two solutions). 
Independently of the terminology, the key reason for estimates (|45|>-(|4"7T> is the Lip- 
schitz property of the operator F. It is remarkable that the large time asymptotics of 
u(x,t), satisfying the problem d28l with such F, can be explicitly expressed through 
spectral characteristics of the linear operator L. 

In order to study the large time asymptotics of u(x,t) in more detail we distin- 
guish two different kinds of asymptotic behavior: 

1) convergence to stationary solutions, 

2) convergence to self-similar solutions provided the condition (c), of the main 
properties on F, is satisfied. 

The case 1) is relatively simple. Any stationary solution u(x) of the problem 
satisfies the equation 



r(u)=u, «eC(R+), ||m|| <1. 



(48) 



If the stationary solution u(x) does exists (note, for example, that F(0) =0 and 
F(l) = 1 for F given in Eqs. (l24l i) then the large time asymptotics of some classes 
of initial data uq(x) in d28l can be studied directly on the basis of Lemma [3] It is 
enough to assume that \uq(x) — u(x)\ satisfies ( PPfl i with p such that A (p) < 1. Then 
u(x,t) — > u{x) as t — > 00, for any x > 0. 

This simple consideration, however, does not answer at least two questions: 
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A) What happens with u(x,t) if the inequality ( l44b for \uq(x) — u(x)\ is satisfied 
with such p that X (p) > 1? 

B) What happens with u(x,t) for large x (note that the estimate (05]) becomes trivial 

if X — > oa). 

In order to address these questions we consider a special class of solutions of 
Eq. (|28K the so-called self-similar solutions. Indeed the property c) of F shows that 
Eq. ( f28l > admits a class of formal solutions u s (x,t) = w(xe^*') with some real jj, t . It 
is convenient for our goals to use a terminology that slightly differs from the usual 
one. 

Definition 3. The function w(x) is called a self-similar solution associated with the 
initial value problem d28l l if it satisfies the problem 

H*9w + w = r(w) , H|<1, (49) 

in the notation of Eqs. j40] |. ( l24b . 

The convergence of solutions u(x,t) of the initial value problem (l28l l to a station- 
ary solution u(x) can be considered as a special case of the self-similar asymptotics 
with jU* = 0. 

Under the assumption that self-similar solutions exists (the existence is proved 
in the next section), we state the fundamental result on the convergence of solutions 
u(x,t) of the initial value problem d28T l to self-similar ones (sometimes called in the 
literature self-similar stability). 

Lemma 4. We assume that 

i) for some ii t £ E, there exists a classical (continuously differentiable if /x* =/= 0) 
solution w(x) ofEq. (|49l such that \\w\\ < 1; 

ii) the initial data u(x,0) = uq in the problem ( |28t satisfies 

uq = w + 0(x p ) , ||w || < 1 ,for p>0 such that ji(p) < ^i*, (50) 

where fl(p) defined in ( 143b is the spectral function associated to the operator L. 
Then 

\u{xe-^ t ,t)-w{x)\ = 0{x p )e- pt ^-^ (51) 

and therefore 

limw(xe^* r ,f) = w(x) , x > . (52) 

Remark 4. Lemma |4] shows how to find a domain of attraction of any self-similar 
solution provided the self-similar solution is itself known. It is remarkable that 
the domain of attraction can be expressed in terms of just the spectral function 
fl(p), p > 0, defined in d43l . associated with the linear operator L for which the 
operator F satisfies the L-Lipschitz condition. 

Generally speaking, the equality d52b can be also fulfilled for some other values 
of p with jli(p) > /I* in Eq. d50l l. but, at least, it always holds if ju(p) < /!*■ 
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We shall need some properties of the spectral function p~(p). Having in mind 
further applications, we formulate these properties in terms of the operator F given 
in Eqs. (I24K though they depend only on K{a) in Eqs. (f42b 

Lemma 5. The spectral function has the following properties: 

i) It is positive and unbounded as p — » + , with asymptotic behavior given by 

M (p)«*MzI, p ^o, (53) 
P 

where, for F from d24b 

„oo A 7 iV 

A(0)= / aaZ(a) = £ <vt > 1 , £ a„ = 1 , a„ > , < 54 ) 

J° n=l n=l 

ana? therefore A (0) = 1 if and only if the operator F (124b ;s linear (N = 1 ); 

ii) /« f/ie case o/a multi-linear F operator, there is not more than one point < ;:>o < °°, 
where the spectral function jJ.(p) achieves its minimum, that is, jU'(po) = 7&{po) 

0, w/f/z ll(po) < ll(p) for any p > 0, provided N > 2 and > 0. 



Remark 5. From now on, we shall always assume below that the operator F from (t24b 
is multi-linear. Otherwise it is easy to see that the problem d49b has no solutions (the 
condition ||w|| < 1 is important!) except for the trivial ones w = 0, 1. 

The following corollaries are readily obtained from lemma[5]part ii) and its proof. 

Corollary 3. For the case of a non-linear F operator, i.e. N > 2, the spectral func- 
tion )l(p) is always monotone decreasing in the interval (0,po)> an d /J-(p) > ^(po) 
for < p < pq. This implies that there exists a unique inverse function p(jtt) : 
(jU(po))+°°) ~~ * (OjPo)> monotone decreasing in its domain of definition. 

Corollary 4. There are precisely four different kinds of qualitative behavior of ju(p) 
shown on Fig.l. 

Proof. There are two options: jJ,(p) is monotone decreasing function (Fig.l (a)) or 
H(p) has a minimum at p = po (Fig.l (b,c,d)). In case Fig.l (a) jj,(p) > for all 
p > since ju(p) > l/p. The asymptotics of X(p) (l42l is clear: 

(1) A(p) >Ao„e]R + if / daK(a)=0; (55) 



(2) > ooif / daK(a)>0. (56) 

In the case (1) when ~ > 00 as P ~ > 0, two possible pictures (with and with- 
out minimum) are shown on Fig.l (b) and Fig.l (a) respectively. In case (2), from 
Eq. d42l) it is clear that X(p) grows exponentially for large p, therefore ji(p) — > °° 
as p — ► oo. Then the minimum always exists and we can distinguish two cases: 
H(po) < (Fig.l (d) and /x(p ) > (Fig.l (c)) 
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Fig. 1 Possible profiles of the spectral function jit (p) 

We note that, for Maxwell models (A), (B), (C) of Boltzmann equation (Sec- 
tions 2 and 3), only cases (a) and (b) of Fig.l can be possible (actually this is the 
case (b)) since the condition (l55l l holds. Fig.l gives a clear graphic representation of 
the domains of attraction of self-similar solutions (Lemma|Ui: it is sufficient to draw 
the line = jU* = constant, and to consider a p such that the graph of ji(p) lies 
below this line. 

Therefore, the following corollary follows directly from the properties of the 
spectral function )J.(p), as characterized by the behaviors in Fig.l, where we assume 
that jii(po) = for po = °o, for the case shown on Fig.l (a). 

Corollary 5. Any self-similar solution u s (x,t) — w{xe >i *') with /J-(po) < fJ.* < °° has 
a non-empty domain of attraction, where po is the unique (minimum) critical point 
of the spectral function fi(p). 

Proof. We use Lemma|4]part ii) on any initial state uq = w + 0(x p ) with p > such 
that jU(po) < M(p) < A 1 *- In particular, Eqs. (f5TJ and d52l show that the domain of 
attraction of w(xe^*') contains any solution to the initial value problem d28| ) with the 
initial state as above. 

It is clear that the inequalities of the kind u\ — U2 = 0(x p ) for any p > such that 
H(p) < /x*, for any fixed /I* > fl(po) play an important role. We can use specific 
properties of jJ.(p) in order to express such inequalities in more convenient form. 
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Lemma 6. For any given /x* G (/i(/?o),°°) and Uip,(x) such that \\u\^\\ < °°, the 
following two statements are equivalent: 

i) there exists p > such that 

u\-u 2 = 0(x p ) , with jx{p) < . (57) 

ii) There exists £ > such that 

mi -u 2 = 0(jc p(/i * )+e ) , w/f/i p( J u*) < po , (58) 
where p(/x) is f/ie inverse to n(p) function, as defined in Corollary|3] 

Finally, to conclude this section, we show a general property of the initial value 
problem d28l for any non-linear F operator satisfying conditions a) and b) given 
in Eqs. d38l ) and (f39l > respectively. This property gives the control to the point-wise 
difference of any two rescaled solutions to ( l28l in the unit sphere of B , whose initial 
states differ by 0(x p ). It is formulated as follows. 

Lemma 7. Consider the problem d28l l. where F satisfies the conditions (a) a«c/ (b). 
Lef Mi 2(^,0 fwo solutions satisfying the initial conditions u\ ,2(^,0) = Mq' 2 (x) 
5mc/z f/zaf 

||mJ' 2 ||<1, u\-ul = 0{xP), p>0. (59) 

f/zen, for any real /i*, 

4 M ,(x,f) = Ul (xe-W,t) -u 2 (xe-^',t) = 0(x p )e- p '^-^ (60) 
and therefore 

]imA fU (x,t)=0 , x>0, (61) 

for any jU* > ju(p). 

Remark 6. There is an important point to understand here. Lemmas [3] and @] hold 
for any operator F that satisfies just the two properties (a) and (b) stated in d38l 
and (|39l ). It says that, in some sense, a distance between any two solutions with 
initial conditions satisfying Eqs. d59l ) tends to zero as t — > °°. Such terminology and 
corresponding distances were introduced for specific forms of Maxwell-Boltzmann 
models in [4, 19]. It should be pointed out, however, that this contraction property 
may not say much about large time asymptotics of u(x,t), unless the corresponding 
self-similar solutions are known, for which the operator F must be invariant under 
dilations (so it satisfies also property (c) as well, as stated in d40b ). In such case one 
can use estimate (|6TT > to deduce the convergence in the form (l52l l. d53l l. 

Therefore one must study the problem of existence of self-similar solutions, 
which is considered in the next Section. 
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5.3 Existence of self-similar solutions 

We develop now a criteria for existence, uniqueness and self-similar asymptotics 
to the problem Eq. ( f49b for any operator F that satisfies conditions (a), (b) and (c) 
from Section 6, with the corresponding spectral function fi(p) defined in d43b . 

Theorem[3]below shows the criteria for existence and uniqueness of self-similar 
solutions for any operator F that satisfies just conditions (a) and (b). Then Theo- 
rem|4]follows, showing a general criteria to self-similar asymptotics for the problem 
d28l l for any operator F that satisfying conditions (a) (b) and (c). 

We consider Eq. d49b written in the form 

Lixw'{x)+w(x)=g(x) , g=r(w) , fieR, (62) 

and, assuming that ||w|| < °°, transform this equation to the integral form. It is easy 
to verify that the resulting integral equation reads 

w(x) = f dxg(xx^) (63) 
Jo 

By means of an iteration scheme, the following result can be proved. 

Theorem 3. Consider Eq. ( 162b with arbitrary ftgl and the operator F satisfying 
the conditions (a) and (b) from Section 6. Assume that there exists a continuous 
function wq(x), x > 0, such that 

i) || wo II < 1 and 

ii) 

/ drgo(x^)=w (x) + 0(xP) , g =r(w ), (64) 
Jo 

with some p > satisfying the inequality 

H(p) = -[[ daK(a)a p — 1] < u . (65) 
p Jo 

Then there exists a classical solution w{x) ofEq. ( |62t . The solution is unique in the 
class of continuous functions satisfying conditions 

\\w\\<l, w(x)=wo(x)+0( X r>) , (66) 

with any p\ such that )J.{pi) < pi- 

Combining Lemma [T] and Lemma |H the following general statement related to 
the self-similar asymptotics for the problem d28l i is obtained. 

Theorem 4. Let u(x,t) be a solution of the problem (128b with \\uq\\ < 1 and F satis- 
fying the conditions (a), (b), (c)from Section 6. Let fl(p) denote the spectral function 
( 1651 ) having its minimum (infimum) at p = pq (see Fig.l), the case po = °° is also 
included. We assume that there exists p 6 (0,/?o) and < £ < Pq— p such that 



20 A.V. Bobylev, C. Cercignani and I.M. Gamba 

[ dT 8o (xT^)=uo(x)+0(x p+e ), go=r(u ) 7 e>0. (67) 
Jo 

Then 

i) there exists a unique solution w(x) of the equation d62t with fl = n(p) such that 

HI < 1 > w ( x ) = uq{x)+0{ X p +£ ) , (68) 

ii) 

lim u(xe^ p) \t ) = w(x) , x > , (69) 
where the convergence is uniform on any bounded interval in R + and 

it(x«-"W,t) - = O^ 6 ^^') , (70) 
withp(p,e) = (p + £)(n(p)-n(p + e))>0. 



Hence, a general criterion ( I68l l is obtained for the self-similar asymptotics of 
m(x,?) with a given initial condition mo(^)- The criterion can be applied to the prob- 
lem d28l l with any operator F satisfying conditions (a), (b), (c) from Section 6. The 
specific class d24l i of operators F is studied in Section 8. We shall see below that the 
condition d68l ) can be essentially simplified for such operators. 



5.4 Properties of self-similar solutions 

We now apply the general theory (in particular, Theorem |4]i to the particular case 
of the multi-linear operators F considered in Section 4, where their correspond- 
ing spectral function ji(p) satisfies ( l65l l. d34l i whose behavior corresponds to Fig.l. 
We also show that po = min p> o ji{p) > 1 is a necessary condition for self-similar 
asymptotics. 

In addition, Theorem [5] establish sufficient conditions for which self-similar so- 
lutions of problem d62l l will lead to well defined self-similar solutions (distribution 
functions) of the original problem after taking the inverse Fourier transform. 

We consider the integral equation ( |63l l written as 

w = r^(w)= [ dt g (xt^) , g = r( w ) , Mem. (71) 
Jo 

The following two properties of w(x) that are independent of the specific form d24l i 
ofF. 

Lemma 8. 

i. If there exist a closed subset U' C U of the unit ball U in B, such that 1^(1/') C U' 
for any ju £ R, and for some function wq £ U' the conditions of Theorem \3\ are 
satisfied, then w £ U'. 
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ii. If the conditions of Theorem\3\for F are satisfied and, in addition, F(l) = 1, 

then the solution w* = 1 ofEq. ( 1711 ) is unique in the class of functions w(x) satisfying 
the condition 

w(x) = l + 0( X P), M(/>)<M- (72) 

We observe that the statement [ii] can be interpreted as a necessary condition for 
existence of non-trivial (w ^ const.) solutions of Eq. ( fTTT i: if there exists a non-trivial 
solution w(x) of Eq. ( fTTT i. where F(l) = 1, such that 

||w|| = l, w=\+0{x p ), p>0, then fi<fi{p). (73) 

We recall that ju(p) satisfies the inequality > n(po) (Fig.l). 
If P > (provided ;;>o < °°) in Eqs. ( 1731 ). then there are no non-trivial solutions 
withjl > n(p ). 

On the other hand, possible solutions with jU < jU(po) (even if they exist) are 
irrelevant for the problem d28l since they have an empty domain of attraction 
(Lemma|4|i. 

Therefore we always assume below that fi> fi (po) and, consequently, p G (0, po) 
in Eq. {73). 

Let us consider now the specific class (|24T> — (f25l of operators F, with functions 
u(x) satisfying the condition m(0) = 1. Then, u(Q,t) = 1 for the solution u(x,t) of 
the problem d281 l. 

In addition, the operators d24"l) are invariant under dilation transformations d40b 
(property (c), Section 6). Therefore, the problem ( 128] ) with the initial condition uq(x) 
satisfying 

a(0) = l, ||«o|| = l; Mo (jc) = 1 -JSjc^H , (74) 

can be always reduced to the case /3 = 1 by transformation x' = xj3 1//p . 

Moreover, the whole class of operators (l24T i with different kernels A„ (a i, ... ,a n ), 
n = 1,2, . . ., is invariant under transformations x = x p , p > 0. The result of such 
transformation of F is another operator F of the same class ( f24b with kernels 

A„(ai,...,a n ). 

Therefore we fix the initial condition d74T i with j3 = 1 and transform the function 
( |74T ) and the equation ( 1281 to new variables x — x p . Then we omit the tildes and re- 
duce the problem ( fZ8b . with initial condition d74l i to the case j3 = 1, p = 1, We study 
this case in detail and formulate afterward the results in terms of initial variables. 

Next, we assume a bit more about the asymptotics of the initial data uq(x) for 
small x, namely 

||«o|| = l, uo{x) = \-x + 0{x l+e ) , x^O, (75) 

with some e > 0. 

Then, our goal now is to apply the general theory (in particular, Theorem |4] and 
criterion d67] i) to this particular case. We assume that the spectral function )J.(p) 
given by d65l l. d34b , corresponds to one of the four cases shown on Fig. 1 with po> 1. 
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Let us take a typical function uq = e~ x satisfying (T73T > and apply the criterion 
d67K from Theorem [4] or, equivalently, look for such p > that d67| i is satisfied. 
That is, find possibles values of p > such that 

WO-^ v = 0(^ +£ ), (76) 

in the notation of Eq. dTTT ). 

It is important to observe that now the spectral function fi (p) is closely connected 
with the operator F (see Eqs. d65l ) and d34li), since this was not assumed in the 
general theory of Sections 4-7. This connection leads to much more specific results, 
than, for example, the general Theorems [3] |4] 

The properties of self-similar solutions to problem d62b for po > 1, and conse- 
quently, for 

Hip) >M(po) > >-l, (77) 

Po 

can be obtained from the structure of r^(e~ x ) for any jU > — 1 using its explicit 
formula (Eqs. (f65t and d34ii). In particular, for 



N 
I 

«=1 • /M '+ u ifc=l 



r ^( e T ) = H a " j^ n dai...da n A n {ai,...,a n )I ii x£a A . 

where 

7 M (y) = f L dte-y" , ;U 6i, y>0, V a„ = 1 



,1 AT 

1 

"l 

the following two statements can be proven. 



(78) 



(79) 



Lemma 9. The condition ( 1761 ) is fulfilled if and only if p < 1, anof therefore fl(p) > 
ju(l) whenever po > 1 wif/z /j(/?o) = rmn p> o 

Theorem 5. 77ie limiting function w(x) constructed by the iteration scheme to prove 
existence in the solution w in Theorem\3\satisfies Eq. ( 1711 ) with /i = jlt(l), where F is 
given in Egs. ( 1241 ), jtt(p) is defined in Eqs. ( 165b , ( 1341 ). 77ien, the following conditions 
are fulfilled for w(x): 

1. 7f satisfies 

<< < 1 , withw(0) = 1 and vv'(0) = -1 , (80) 
w'{x)<0, \w'(x)\ < 1 , ant/ w(jc) = e~ x + 0{^) , (81) 



wif/i 

7r(jli) = <J 2 e foranye>0 ifn = ±, 

p i/--l<M<-i 



2. Further 
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e~ x < w(x) < 1 , lim w(x) = , and (82) 
3. there exists a generalized non-negative function R(x), T > 0, such that 

w(x)= [ d%R(x)e- xx , f dxR(x) = f dxR(x)x = l. (83) 
Jo Jo Jo 

The integral representation d83l is important for properties of corresponding dis- 
tribution functions satisfying Boltzmann-type equations. Now it is easy to return to 
initial variables with uq given in Eq. d74) and to describe the complete picture of the 
self-similar relaxation for the problem d28) . 



6 Main results for Maxwell models with multiple interactions 
6.1 Self -similar asymptotics 

We apply now the results of Section 6 to the specific case when the Cauchy problem 
d28l ) with a fixed operator F (124") corresponds to the Fourier transform problem for 
Maxwell models with multiple interactions. In particular we study the time evolu- 
tion of uq(x) satisfying the conditions 

IKH = 1 ; u = l-x p + 0(x p+£ ) , x -> , (84) 

with some positive p and e. Then, from Theorems [3] and HI there exists a unique 
classical solution u(x,t) of the problem d28l ). d84l ) such that, for all t > 0, 

IK-,f)ll = 1 ; u{x,t) = l + 0{x p ) , x->0 . (85) 

First, consider the linearized operator L given in Eqs. d33l)-(l34l) and construct the 
spectral function ji (p) given in Eq. d65l l which will be of one of four kinds described 
qualitatively on Fig. 1 . 

Second, find the value po > that corresponds to minimum (infimum) of fi(p). 
Note that po = °° just for the case described on Fig.l (a), otherwise < po < °°. 
Compare po with the value p from Eqs. (f84T >. If p < po then the problem d28l i. ([84-b 
has a self-similar asymptotics (see below). 

In particular, two different cases are possible: ( 1 ) p > po provided po < °°; (2) < 
p < po- In the first case a behavior of u(x,t) for large t may depend strictly on initial 
conditions. 

Depending on how p compares with po, we can obtain We again use Lemma [7] 
with u\ = u and «i = u s = xjfixe^P^') and obtain for the solution u(x,t) of the prob- 
lem (EH), (El): 
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[l if M>M(P) 

lim u(xe-V,t) = { w(x) if fi = n{p) (86) 

[0 if n{p)>y,>li{p + 8) , 

with sufficiently small 8 > 0. 

We see that = where w(x) has all properties described in Theorem|5] 

The equalities d86t explain the exact meaning of the approximate identity, 



u(x,t) w y/ixe^') , t -> oc , xe^'' = const. , (87) 

that we call self-similar asymptotics. 

In particular, the following statement holds. 

Proposition 1. The solution u(x,t) of the problem ( 128b , (1841 1. with r given in 
Eqs. ( 1241 ), satisfies either one of the following limiting identities: 

(1) 

lim M (xe^',f) = 1 , x>0, (88) 



for any jl > /J.{po). if P > Pofor the initial data 
(2) Eqs. (l86l l provided < p < po. 

The convergence in Eqs. ( 1881 , d86t w uniform on any bounded interval < x < R, 
and 

u(xe^',t) - W {x) = 0(x? +£ )e-P^' , j5(p, e) = (p + e) (/x(p) -M(p + e) ), 
/or < p < po and < E < po — p. 

Remark 7. There is a connection between self-similar asymptotics and non-linear 
wave propagation. It is easy to see that self-similar asymptotics becomes more trans- 
parent in logaritheoremic variables 

y = Inx , u(x,t) = u(y,t) , ys(x,t) = fy(y,t) 

Thus, Eq. j87l becomes 

u(y,t) w f(y + n(p)t) , t-K*>, y + n(p)t = const. , (89) 

hence, the self-similar solutions are simply nonlinear waves (note that = 1, 

= 0) propagating with constant velocities c p = —fi (p) to the right if c p > 
or to the left if c p < 0. If c p > then the value u(— °°,t) — 1 is transported to any 
given point y g R when f — > °°. If c p < then the profile of the wave looks more 
naturally for the functions u = 1 — u, \fr = 1 — yf. 

We conclude that Eq. d28b can be considered in some sense as the equation for 
nonlinear waves. The self-similar asymptotics d89b means a formation of the trav- 
eling wave with a universal profile for a broad class of initial conditions. This is a 
purely non-linear phenomenon, it is easy to see that such asymptotics cannot occur 
in the particular case (N = 1 in Eqs. ( l24l ) of the linear operator F. 
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6.2 Distribution functions, moments and power-like tails 



We have described above the general picture of behavior of the solutions u(x,t) to 
the problem (l28l l. (|84| |. On the other hand, Eq. ( T28l > (in particular, its special case 
O) was obtained as the Fourier transform of the kinetic equation. Therefore we 
need to study in more detail the corresponding distribution functions. 

Set uo(x) in the problem (l28l l to be an isotropic characteristic function of a prob- 
ability measure in R d , i.e., 

«o(jc) = &\fo] = [ dvf (\v\)e- ik v , k e R d , x = \k\ 2 , (90) 

where fo is a generalized positive function normalized such that uq(Q) = 1 (distri- 
bution function). Let U be a closed unit ball in the B = C(M+) as defined in (f30b . 

Then, we can apply all results of Section 5, and conclude that there exists a 
distribution function f(v,t), v £ K d , satisfying Eq. (|T), such that 

U(x,t)=&\f(;t)], X=\k\\ (91) 

for any t > 0, and a similar conclusion can be obtain if we assume the Laplace 
(instead of Fourier) transform in Eqs. (1901 1. 

Then there exists a distribution function f(v,t ), v > 0, such that 

u(x,t) =W(; t)] = [ dvf(v,t)e- xv , k(0,0 = 1, x>0, t>0, (92) 
Jo 

where u(x,t) is the solution of the problem (l28T > constructed in Theorem |2] and 
Lemma|2] 

The approximate equation (f87T > in terms of distribution functions (PTT i reads 

f(\v\,t)~e-^(p)>F p (\v\e-^(p)t) ; f^oo, IvIe-sfW' = const. , (93) 
where F p (|v|) is a distribution function such that 

yr P (x)=&[F p ], x=\k\ 2 , (94) 

with y/p given by 

it^x,*) = Vf(*"W) (95) 

(the notation is used in order to stress that y/ defined in d95l l. depends on p). 
The factor 1 /2 in Eqs. d93l is due to the notation x = \k\ 2 . Similarly, for the Laplace 
transform, we obtain 

f(v,t) ~ e-^'cppive-^') , t -> oo , ve -M(rt' = const i (96 ) 

where 

yr p (x) =l[* p ] . (97) 
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The positivity and some other properties of F p (|v|) follow from the fact that 
typ(x) = w p (x p ), where w p (x) satisfies Theorem|5] Hence 

W P {x) = [°" dxR p {x)e- txP , f°°dxR p (x) = [°°dxR p {x)x = 1 , (98) 
Jo Jo Jo 

where R p (x), X > 0, is a non-negative generalized function (of course, both \j/ p and 
R p depend on p). 

In particular we can conclude that the self-similar asymptotics ((93) for any initial 
data /o > occurs if po > 1, otherwise it occurs for p £ (0,/?o) C (0, 1 ). Therefore, 
for any spectral function jJ.(p) (Fig.l), the approximate equality (|93l holds for suf- 
ficiently small < p < 1 . In addition, 

m 2 = / c/v/o(|v|)|v| 2 < °o if p = 1 

and 1112 = °° if p < 1. Similar conclusions can be made for the Laplace transforms. 

The positivity of F( |v|) in Eqs. d94ll-(l97]i follows from the integral representation 
d98t with p < 1, since it is well-known that 

^-i( e -l*l 2 ") >0 , L- 1 (e-^)>0 

for any < p < 1 (the so-called infinitely divisible distributions [18]). Thus, 
Eqs. d98l l explains the connection of the self-similar solutions of generalized Maxwell 
models with infinitely divisible distributions. 

Using standard formulas for the inverse Fourier (Laplace) transforms, denote by 
(d = 1,2, . . . is fixed) 

(99) 

1 ra-\-ioa 

N P (y) = T-- ■ dxe ~ xP+xv . 0<p<l, 

2,7^1 J ci — foo 

we obtain the self-similar solutions (distribution functions) are given in Eqs. d96l ), 
d98l ) (right hand sides), by 

F p (|v|) = J dxR p {x)x-i'M p {\v\x-^) , 

<P p (v) = J dTR p {T)T^m p {vx^p) , V>0,0</?<1. 

Note that Mi(|v|) is the standard Maxwellian in R d . The functions N p (v) d99l are 
studied in detail in the literature [18,20]. Thus, for given < p < 1, the kernel 
R(x), X > 0, is the only unknown function that is needed to describe the distribution 
functions F(|v|) and 4>(v). See [10] for a study of R(x) in more detail. 
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Now, from (|34] >. ( l42l . and recalling jj. = fx (1 ), we can show the moments equation 
can be written in the form 

N 

(sn(l) - A (j) + l)m, = ^ . (101) 

n=2 

where 

/„(*) = /c/fli ...da„A(a\,...,a„) dxi...,dx„ g^\aiXi,...,a n x n )Y\R(Xj) 



^ ) (yi ) ...,y») = ( 1%) -£)£, » = i,2,. 

A=l / fc=l 



(102) 



and, due to the properties of R(x), one gets = for any s > and mo = wi = 1. 

Our aim is to study the moments m s defined in Eq.(|92l, for s > 1, on the basis of 
Eq. d 1 1 b - The approach is similar to the one used in [15] for a simplified version of 
Eq. ( 11011 ) with N = 2. The main results are formulated below in terms of the spectral 
function ji(p) (see Fig.l) under assumption that po > 1. 

Proposition 2. 

[i] If the equation fl(s) = /i(l) has the only solution s = 1, then m s < °°for any s > 0. 
[ii\If this equation has two solutions s = 1 and s = 5* > 1, then m s < °°for s < s t . and 

m s = °°for s > i*. 
[iiiKs, < °° only ifl„(s. t ) = in Eq. ( 1101b for all n = 2, ... ,N. 

Now we can draw some conclusions concerning the moments of the distribution 
functions dlOOb as follows. Denote by 



m s (%)= dv%(v)v s , m s {R p )= drR p (T)T s , 
Jo Jo 

m 2s (F p ) = f dvF p (\v\)\v\ 2s , s > , < p < 1 , 

and use similar notations for N p (v) and M p (|v|) in Eqs. (1100b . Then, by formal 
integration of Eqs. (1 100b . we obtain 

m s {<P p ) = m s (Np)m s i p (Rp) 

m2s(Fp)=m2s(M p )m s / p (Rp) , 



where M p and N p are given in Eqs. ( f99b and we show that finite only for s < p < 1. 
In the remaining case p = 1 all moments of functions 



28 



A.V. Bobylev, C. Cercignani and I.M. Gamba 



Mi(|v|) = (4*) 



exp 




,d . 



Ni(v) = S(v-l) 



rel 



■+ > 



are finite. Therefore, everything depends on moments of R\ in with p = 1, so it only 
needs to apply Proposition^ 

In particular, the following statement holds for the moments of the distribution 
functions d93l , d96l l. 

Proposition 3. 

[i] I/O < p < 1, then m2 s (F p ) and m s {<t> p ) are finite if and only ifO < s < p. 
[ii\If p = 1, then Proposition\2\holds for m s = ni2 s {F\) and for m s = m s (<Pi). 

Remark 8. Proposition[3]can be interpreted in other words: the distribution functions 
F p (|v|) and 4> p (v), < p < 1, can have finite moments of all orders in the only case 
when two conditions are satisfied 

(1) p = l,and 

(2) the equation fi(s) = (see Fig.l) has the unique solution s = 1. 

In all other cases, the maximal order s of finite moments ni2 S (F p ) and m s (<P p ) is 



This fact means that the distribution functions F p and <P p have power-like tails. 



6.3 Applications to the conservative or dissipative Boltzmann 
equation 

We recall the three specific Maxwell models (A), (B), (C) of the Boltzmann equation 
from section 2. Our goal in this section is to study isotropic solutions /(|v|,f),v G 
K d , of Eqs. (O, ©, and © respectively. All three cases are considered below from 
a unified point of view. First we perform the Fourier transform and denote 



It was already said at the beginning of section 4 that u(x,t) satisfies (in all three 
cases) Eq. (l23l . where = 2 and all notations are given in Eqs.d27l>. ([Hi t. ([TUt-lETTi. 
Hence, all results of our general theory are applicable to these specific models. In 
all three cases (A), (B), (C) we assume that the initial distribution function 



bounded. 




(103) 





(104) 



and the corresponding characteristic function 
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u(0,t)=u o (x)=^\f (\v\)], x=\k\\ (105) 

are given. Moreover, let uq(x) be such that 

uq(x) = 1 - ax p + 0(x p+e ), x^O, 0<p<l, (106) 

with some a > and e > 0. We distinguish below the initial data with finite energy 
(second moment) 

E = [ rfv|v| 2 /o(|v|)<°° (107) 



implies p = 1 in Eqs. ( 1106b and the in-data with infinite energy Eq = °°. If p < 1 in 
Eqs. < PT05t then 

mf= [ dvf (\v\)\v\ 2 «<~ (108) 

only for 4 < p < 1 (see [18,20]). 

Also, the case p > 1 in Eqs. ( 11061 ) is not possible for /o(|v|) > 0. 

In addition, note that the coefficient a > in Eqs. ( 11061 ) can always be changed 
to a = 1 by the scaling transformation x = a}l p x. Then, without loss of generality, 
we set a = 1 in ( 1106b . 

Since it is known that the operator F(m) in all three cases belongs to the class 
((24b . we can apply Theorem [5] and state that self-similar solutions of Eq. d23l l are 
given by 

u s {x,t) = f (xe^'), f (jc) = w{x p ) , (109) 

where is given in theorem|5]and < p < po (the spectral function fi (p), defined 
in (|43t . and its critical point p$ depends on the specific model.) 

According to Sections 9-10, we just need to find the spectral function fx(p). In 
order to do this we first define the linearized operator L = -T'(l) for r(u) given in 
Eqs. (1241) . (I27t . One should be careful at this point since A2 (a 1,122) in Eqs. (l27b is not 
symmetric and therefore Eqs. d32t cannot be used. A straight-forward computation 
leads to 



Lu(x) = / dsG(s)(u(a(s)x)+u(b(s)x))+ dsH(s)u(c(s)x) , (110) 
Jo Jo 

in the notation ( fT8l - d2U . Then, the eigenvalue X(p) is given by 
Lx p = X(p)x p which implies 

k{p)= f l dsG(s){(a(s)) p + (b(s)) p }+ [ l dsH(s){c(s)) p , 
JO Jo 

and the spectral function (l43l l reads 

M(P) = ^. dl2) 
Note that the normalization (f22b is assumed. 
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At that point we consider the three models (A), (B), (C) separately and apply 
Eqs. ( Ill U and (II 12b to each case. 

(A) Elastic Boltzmann Equation © in R d ,d > 2. By using Eqs. (fT8T >. (fl9l and d22ji 
we obtain 

A(p) = [ 1 dsG(s)(sP + (l-s)P), G{s)=A d g(l-2s)[s(l (113) 
Jo 

where the normalization constant A,./ is such that Eq. (l22l is satisfied with H = 0. 
Then 



1 /-i 

M(/0 = - / *Gf(* i P + (l-jj'-l) ,p>0. 
p Jo 



r 

p. 

It is easy to verify that p/i(p) — > 1 as p — > 0, jlt(p) — > as p — > «>, and 



(114) 



ju(p) > if p < 1 ; ju(p) < if p > 1 ; 

rl (115) 
M(1)=0, /i(2) = /i(3) = - / ^G(i)i(l-*). 



Hence, /I (p) in this case is similar to the function shown on Fig. 1 b) with 2 < po < 3 
and such that =0. Then the self-similar asymptotics hold all < p < 1. 
(B) Elastic Boltzmann Equation in the presence of a thermostat © in R d ,d > 2. 
We consider just the case of a cold thermostat with T = in Eq. (fl4l . since the 
general case T > can be considered after that with the help of (fTTT i. Again, by 
using Eqs. (fT8l , (fT~9b , and (l22l we obtain 



A(p)= / dsG(s)(sP + (l-sy) + G ( dsG( S )(l-—^y, 
Jo Jo l+m) z 



G(s) = T ^A d g(l-2 S )ls(l-s)} £ 



(116) 



with the same constant Aj as in Eq. (II 13b - Then 

M(p) = -/ 1 <faG(*)(y + (l- J )'-0(l-W-(l + 0)) ] 
p Jo 



4m 

(iT^) 2 



(117) 



and therefore, as in the previous case A), p /I (p) — * 1 as p — > 0, ji (p) — > as p — > oo, 
with 

u(l) = -0fl [ 1 dsG(s)s. (118) 
Jo 

which again verifies that ji (p) is of the same kind as in the elastic case A) and shown 
on Fig.l (b). A position of the critical point po such that /i'(po) = (see Fig.l (b)) 
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depends on 9. It is important to distinguish two cases: 1) po > 1 and 2) po < 1. In 
case 1) any non-negative initial data ( 1104b has the self-similar asymptotics. In case 
2) such asymptotics holds just for in-data with infinity energy satisfying Eqs. ( 1106b 
with some p < po < 1 . A simple criterion to separate the two cases follows directly 
from Fig.l (b): it is enough to check the sign of jtt'(l). If 

(i'(l)=l'(l)-l(l)|l<0 (119) 

in the notation of Eqs. dl 16b , then po > 1 and the self-similar asymptotics hold for 
any non-negative initial data. 

The inequality dl 19b is equivalent to the following condition on the positive cou- 
pling constant 9 

0<fl<a> = _/>^) (^+(1-^(1-^)), (120) 

/ ( >G(*)(j3*+(l-j3*)log(l -/?*)) 

The right hand side of this inequality is positive and independent on the normaliza- 
tion of G(s), therefore it does not depend on 9 (see Eq. dl 17b . We note that a new 
class of exact self-similar solutions to Eq. (O with finite energy was recently found 
in [12] for j3 = 1, 9 = 4/3 and G(s) = const. A simple calculation of the integrals 
in ( 1120b shows that 0* = 2 in that case, therefore the criterion dl 19b is fulfilled for 
the exact solutions from [12] and they are asymptotic for a wide class of initial data 
with finite energy. Similar conclusions can be made in the same way about exact 
positive self-similar solutions with infinite energy constructed in [12]. Note that the 
inequality ( 11 19b shows the non-linear character of the self-similar asymptotics: it 
holds unless the linear term in Eq. (@]i is 'too large'. 

(C) Inelastic Boltzmann Equation (O in W 1 . Equations (fJB and d22l lead to 



X(p)= f dsG{s){{asY + {\-bs)P). 
Jo 



where 

« (l+e) 2 u (l+e)(3 



G(s)=C d (l-sy 



4 ' 4 
with such constant Q that Eq. ( 1221 with H = is fulfilled. Hence 

U(p) = - f 1 dsG(s)((asy + (l-bs) p -l), p>0, (121) 
p Jo 

and once more as in the previous two cases, p jj, (p) — > 1 as p — ► 0, fi (p) — > as 
p — > °°, with now 

1-e 2 r l 



M(l) = 7^ / dsG(s)s. 

4 Jo 
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Thus, the same considerations lead to the shape of ju(p) shown in Fig.l (b). The in- 
equality dl 19b with X(p) given in Eqs. (1121b was proved in [11] (see Eqs. (4.26) 
of [11], where the notation is slightly different from ours). Hence, the inelastic 
Boltzmann equation (0 has self-similar asymptotics for any restitution coefficient 
< e < 1 and any non-negative initial data. 

Hence, the spectral function /J.(p) in all three cases above is such that po > 1 
provided the inequality ( 1120b holds for the model (B). 

Therefore, according to our general theory, all 'physical' initial conditions d 1 041 > 
satisfying Eqs. ( 1106b with any < p < 1 lead to self-similar asymptotics. Hence, the 
main properties of the solutions f(v,t) are qualitatively similar for all three models 
(A), (B) and (C), and can be described in one unified statement: Theorem|6]below. 

Before we formulate such general statement, it is worth to clarify one point re- 
lated to a special value < p\ < 1 such that fi(pi) = 0. The reader can see that on 
Fig. 1 (b) that the unique root of this equation exits for all models (A), (B),(C) since 
jll(l) = in the case (A) (energy conservation) , and < in cases (B) and (C) 
(energy dissipation). If p = p\ in Eqs. ( 1106b then the self-similar solution d 1 09b is 
simply a stationary solution of Eq. d23l ). Thus, the time relaxation to the non-trivial 
(u ^ 0, 1) stationary solution is automatically included in Theorem|6]as a particular 
case of self-similar asymptotics. 

Thus we consider simultaneously Eqs. (|2j, ©, with the initial condition 
( 1104b such that Eq. ( 1106b is satisfied with some 0<p<l,£>0 and a = 1. We 
also assume that T = in Eq. (|4]i and the coupling parameter 6 > satisfies the 
condition ( 11201 ). 

In the following Theorem|6l the solution /(|v|,f) is understood in each case as a 
generalized density of probability measure in Mr and the convergence /„ — > / in the 
sense of weak convergence of probability measures. 

Theorem 6. The following two statements hold 

[i] There exists a unique (in the class of probability measures) solution f(\v\,t) to each 
of Eqs. (O, (|4), © satisfying the initial condition ( 1104b . The solution /(|v|,f) has 
self-similar asymptotics in the following sense: 

For any given < p < 1 in Eqs. ( 11061 ) there exits a unique non-negative self-similar 
solution 

f( p \\v\,t)=e-^ t F p (\v\e-^ t ), (122) 

such that 

e^'f(\v\e-^',t) ^_F p (|v|), (123) 

where is given in Eqs. (II 14b . (II 17b , (1121b . respectively, for each of the three 
models. 

\ii\Except for the special case of the Maxwellian 

F x {\v\) =M(|v|) - (4n)- d/2 e-^ (124) 
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for Eq. (|2]i with p = 1 in Eq. ( 1106t (note that ju(l) =0 in f/iis case), the function 
F p {\v\) does not have finite moments of all orders. IfO < p < 1, then 



If p = 1 in the case ofEqs. ©, ©, f/ien m ? < °° onZy /or < q < p t , where p* > 
1 is the unique maximal root of the equation )J.(p*) = Jlt(l), with ll(p) given in 
Eqs. ( 1106b , (112U respectively. 

In addition, we also we obtain the following corollary. 

Corollary 6. Under the same conditions of Theorem® the following two statements 
hold. 

[i] The rate of convergence in Eq. (1123b is characterized in terms of the corresponding 
characteristic functions in Proposition^ 

[ii] The function F p (\v\) admits the integral representation ( 1 1001 ) through infinitely di- 
visible distributions ( 1991 . 

Proof. It is enough to note that all results of section 9 and 10 are valid, in particular, 
for Eqs. ©, ©, ©. 

Finally, we mention that the statement similar to theorem[6l can be easily derived 
from general results of sections 7.1 and 7.2 in the case of 1 — d Maxwell models 
introduced in [15], [2] for applications to economy models (Pareto tails, etc.). The 
only difference is that the 'kinetic' equation can be transformed to its canonical form 
d23l- (l24l i by the Laplace transform as discussed in section 10, and that the spectral 
function ji (p) can have in this case any of the four kind of behaviors shown in Fig. 1 . 
The only remaining problem for any such 1— dimensional models is to study them 
for their specific function ji(p), and then to apply Propositions Q] [2] and [3] 

Thus, the general theory developed in this chapter is applicable to all existing 
multi-dimensional isotropic Maxwell models and to 1 d models as well. 
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